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Abstract 

In grand unified theories based on orbifold constructions in higher dimensions, Higgsino- 
mediated proton decay is absent. However, proton decay mediated by X and Y gauge 
bosons is typically enhanced to levels detectable by current and future experiments. We 
analyse the phenomenology of proton decay induced by the minimal coupling of X, Y 
gauge bosons. In particular, we show that the novel realization of matter in orbifold 
GUTs can lead to unusual final state flavour structure, for example, the dominance of 
the p —>■ K^fi~^ mode. Furthermore, we discuss proton decay induced by higher-derivative 
brane operators, finding potentially observable rates for natural values of the operator 
coefficients. 



1 Introduction 



Grand unified theories (GUTs) provide an elegant explanation of the origin of the three 
Standard Model (SM) gauge interactions and the fermion quantum numbers In their 
minimal supersymmetric extension, GUTs also lead to a remarkably successful prediction 
of sin^ 6w. This supports both the existence of supersymmetry (SUSY), broken near the 
weak scale, and some form of gauge-coupling unification at a scale Mqut — 10^^ GeV. 
The idea of GUTs can also naturally accommodate relations among Yukawa couplings, 
leading most notably to the quite successful prediction of mb/mr- 

Among the less attractive features of GUTs are the comphcated Higgs sector required 
for realistic breaking of the gauge group and the necessity of modifying the unsuccess- 
ful first and second generation analogues of the mb/mj. mass predictions. In addition, 
dimension-5 proton decay operators arising from the exchange of supermassive coloured 
Higgsinos severely constrain minimal SUSY GUT models (see e.g. 0). 

Recently an elegant solution to these problems has been proposed in the context of 
SU(5) and SO(IO) unification. The GUT gauge symmetry is now realized in 
5 or more space-time dimensions and broken to the SM group by compactification on 
an orbifold, utilizing boundary conditions that violate the GUT-symmetry. In the most 
studied case of 5 dimensions both the GUT group and 5d supersymmetry are broken 
by compactification on jiZ^ x Zg), leading to an N=l SUSY model with SM gauge 
group. This construction allows one to avoid some unsatisfactory features of conventional 
GUTs with Higgs breaking, such as doublet-triplet splitting, dimension-5 proton decay, 
and Yukawa unification in the first two generations, while maintaining, at least at leading 
order, the desired gauge coupling unification P,[5|,pTi|. 



Higgsino mediated proton decay is absent because the triplets acquire mass via the 
KK expansion of 5d terms of the form 

£ D ^ SB [ny^Hs + H^V^H^) + h.c. (1) 

Thus, the mass couples the triplet Higgsino to a state in H^, and not in H-^. Unlike 
the H3 and H-^ states, the and if| fields do not couple directly to the quark and 
lepton superfields and the dangerous dimension-5 operators are absent. This absence 
can be viewed as a consequence of a U(l)ij symmetry of the 5d theory. An extension of 
the U(l)/j to matter fields localized on the branes leads to i?-parity, prohibiting baryon 
number violation at dimension 4 as well 

Although there is no dimension-4 or -5 proton decay, we argue in this letter that 
X, Y gauge-boson mediated proton decay is now much more interesting. First, the mass 
scale Mc = 1/R oi the X and Y gauge bosons in the orbifold GUT theories is lower 
than in conventional SUSY GUTs, possibly approaching 10^^ GeV.Q This enhances the 
dimension-6 proton decay processes to a level that may be seen in current and future 
experiments. Second, the novel realization of the matter multiplets in orbifold GUTs 

^The reason for this is that the running of differences of gauge couphngs does not stop at Afc, but 
continues in only a sHghtly modified way |||,||]. So Mc is lower than the unification scale. 
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changes in a striking way the signatures of X, Y mediated proton decay and leads to a 
clear experimental distinction between orbifold and 4d GUT predictions.^ 

We begin in Sect. 2 by deriving the dimension-6 proton decay operators arising from 
the minimal couplings of X, Y gauge bosons. We integrate out the 5d X, Y states in a 
4d superfield formalism. The advantage of this approach is that no ill-defined contact- 
interaction ~ 5(0) (cf. ||12|) appears in the intermediate stages of the calculation. Using 
this formalism we then discuss the effect of higher-derivative brane operators involving 
the X, Y gauge fields. Such interactions, which are not forbidden by the gauge symme- 
tries of the model, can lead to sizeable proton decay rates even in scenarios which are 
otherwise safe due to the location of the matter multiplets. In Sect. 3, we calculate the 
corresponding proton decay rates for the two classes of operators. The obtained results 
depend crucially on the location of matter multiplets on the two branes or in the bulk. 
In particular, we point out that in one of the most attractive realizations of the SU(5) 
orbifold model, the minimal X, Y couplings favour p — »• n'^K^ decay. This represents 
a striking orbifold GUT signature. We also show that, depending on the realization of 
matter, the higher-derivative operators involving X, Y gauge bosons, can be the numer- 
ically dominant source of proton decay. Overall, X, F-mediated proton decay can occur 
at rates observable in the current or next generation of proton decay experiments. Our 
conclusions are given in Sect. 4. 



2 Integrating out the X, Y gauge bosons 

We begin by recalling the basic structure of the Kawamura model , which is based on 
a 5d super Yang-Mills theory on iR^ x S^, where the is parameterized by y G [0, 2ttR). 
The field space is then restricted by imposing the two discrete Z2 symmetries, y —>■ —y 
and y' —y' (with y' = y — ttR/2). The action of the Z2's in field space is specified by 
the two gauge twists P and P' . If the original gauge group is SU(5) and the gauge twists 
are chosen as P = 1 and P' = diag(l, 1, 1, —1, —1), an effective low energy theory with 
SM gauge symmetry results. 

The supersymmetric version of this model can be described in terms of 4d super- 
fields ||13| . A manifestly gauge- invariant form of the non-abelian action is given by |14 



where \^ is a real superfield depending on the additional parameter y = x^, W is the 
corresponding field strength superfield, V5 = c^s + $ is the covariant derivative in 
direction, and $ is an x^-dependent chiral superfield. 

We will be interested in baryon-number-violating processes mediated by X, Y gauge 
bosons in scenarios where SM fermions are localized on the SU(5) brane at y = 0, on 
the SM brane aX y = I = 7ri?/2, or in the bulk. We begin with the simplest case of 

^We thank Y. Nomura for bringing to our attention Ref. |Q where these interesting features of 
orbifold GUT theories were previously noted. In Sect. 3 we will comment on the relations of our respective 
findings. 
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SU(5)-brane localized matter. The coupling to a chiral superfield at y = is described 
by the lagrangian 

= 5{y) I ^e^^v^ = 5{y) f + 2V^J^ + 0{V^)] , (3) 

where =tr(\l/T'^\&) is the matter current. Note that the SU(5) generators {T^} fall 
into the two subsets {T"} and {T°}, the first one being the SM generators and the second 
the basis of the orthogonal complement. 

Given a current J (corresponding to an external field \E'), the equation of motion 
for V is obtained by varying the action of Eqs. (0) and (^. In the $ = gauge and to 
leading order in V it reads (our superspace conventions are those of |[T5|] ) 



- ^D^d'd'^V - dlV + glJb{y) = . (4) 

If is on-shell, Z)^\I' = 0, one finds D^J = 0. It is self-consistent to assume that, in the 
presence of an x'^-independent current (/i = 0, .., 3), we have an independent solution 
V with a superspace dependence given by \^ ~ J. Then the first term in Eq. (^) vanishes 
and we need to solve 

- dIV' + glrSiy) = (5) 

for the V components corresponding to X, Y gauge bosons. Given the boundary con- 
ditions d^V^iy = 0) = V^{y = I) = 0, the solution is = -gl{l - y)J^ for y > 0. 
The function is continuous at ?/ = but has a discontinuous first derivative. Inserting 
this into the original lagrangian, the following effective operator in the low-energy 4d 
effective theory is generated: 

Oi = 'i9ll^-,T.(rf- (6) 

a 

We have checked explicitly that the four-scalar interaction contained in the superfield op- 



erator Oi agrees with the result of , where the calculation was done by first integrating 
out auxiliary fields and then summing the exchanged KK modes in the component formal- 
ism. The advantage of our superfield approach is that no ill-defined contact-interaction 
~ 5(0) appears in the intermediate stages of the calculation^. This could also have been 
achieved by summing KK modes and making use of the cancellations between physical 
and auxiliary field contributions enforced by the unbroken 4d SUSY at each KK level. 

The four-fermion interaction contained in Oi reads 



c a 



(7) 



where ip is the fermion from \E', gl = g^/l is the 4d gauge coupling, and Mc = l/i? is 
the compactification scale which, at the same time, is the mass of the lowest-lying KK 
mode of X,Y gauge bosons. Note that, up to the prefactor 7r^/4, this is precisely the 
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This observation has also been made in the context of a 3d superfield description of 4d SUSY 



theories with boundary |16 
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result expected in a 4d GUT with gauge boson mass M^. The prefactor can be easily 
understood as 2 (vr^/S), where the 2 comes from the normalization of the KK modes due 
to their non-trivial bulk profile and vr^/S = I]^i(2n — accounts for exchanging 
the full KK tower rather than just the lowest-lying mode. 

We now turn to the case where SM fermions are localized on the SM brane 
Although in this case, X, Y mediated proton decay is no longer a generic prediction of 
the theory, it will occur if brane operators of the type 



U = ^Siy') I M/i (Vge^^) ^2 + h.c. 



are included in the action |]T4|. Here and '^2 are SM superfields and the multiplication 



with (VsC^^) is defined by using their standard embedding into SU(5) multiplets. Note 
that, even though the X, Y components of V vanish at the SM brane, their c^s derivatives 
appearing in (Vse^^) are non-zero. The prefactor includes a dimensionless coeffi- 
cient c and the fundamental scale M ^ M,,, which is required for dimensional reasons. 
Parametrically, M can not be larger than ~ l/fi'5, since this is the scale at which the 5d 
gauge theory becomes strongly coupled. 

Following the same line of reasoning as in the case of SU(5)-brane fermions, one 
arrives at the analogue of Eq. (^) , which is 

-55V^-^^'55%)=0, (9) 



with the current J" =tr(\E'iT'^\E'2+h.c.). The solution is given by 

ly 



V\y) = ^rfdy5{y). (10) 



Inserting this into the original lagrangian, the following effective operator in the low- 
energy 4d effective theory is obtained: 

a '-a 

Here, to obtain the final expression, we have assumed the 5 function to be localized on 
a scale 1/M and, accordingly, replaced the y integral by 6M, with h an 0{1) coefficient 
depending on the specific way in which the 5 function is regularized. Thus, the operator 
O2 depends via on the a priori unknown strength of the brane-bulk coupling, Eq. (|]), 
and via h on the UV completion of the theory. This UV sensitivity can also be understood 
by observing that the derivative coupling of Eq. (H) enhances the contribution of higher 
KK modes so that the sum diverges. Parametrically, O2 is suppressed with respect to 
Oi by a factor M/Mc-, which can be as small as ~ 10 (allowing at least some range of 
validity for the effective 5d theory) or as large as ~ 10^ (extending the 5d theory all the 
way to its strong coupling scale). Nevertheless, as we argue in the next section, in models 
where the first generation is not on the SU(5) brane, O2 can be a competitive or even 
the dominant source of proton decay. 

^Following such realizations of matter have also been used in recent TeV scale SU(3) models JiTt. 
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In models where both branes contain SM fermions, the combined effect of the cou- 
phngs of Eqs. (|^) and (§) will give rise to further proton decay operators involving both 
SU(5)- and SM-brane fermions. The calculation is a straightforward combination of the 
two basic cases discussed above. 



3 Proton decay rates 

In SU(5) orbifold GUTs, there are three possible locations for matter: the SU(5) brane, 
the bulk (which is also SU(5) symmetric), and the SM brane. The various possible models 
are characterized by the placement of the Tj (lO's) and Fi (5's) that make up the three 
{i = 1, 2, 3) generations of quarks and leptons. To reproduce the successful SU(5) mt/m^ 
mass relation, T3 and F3 must reside on the SU(5) brane. Furthermore, due to the 
different normalization of 4d and 5d fields, Yukawa interactions involving one or two bulk 
matter fields are effectively suppressed by factors (M/Mc)^/^ and (M/Mc) (recall that 
M is the fundamental UV scale of the theory). Thus, both the relative lightness of the 
2nd and 1st generations and the failure of the SU(5) mass predictions indicate that for 
i = 1,2 either Tj or Fi or both should be in the bulk or on the SM brane. 

First we consider the case of proton decay operators arising from minimal X, Y gauge 
boson interactions. 

Both for matter on the SM brane and in the SU(5)-symmetric bulk, minimally cou- 
pled X, Y gauge bosons do not lead to dimension-6 baryon number violating opera- 
tors [ITT], 113. For matter on the SM brane this is simply because the orbifold boundary 



conditions imply the vanishing of the X, Y gauge boson wavefunctions on the SM brane. 
For bulk matter the situation is slightly more involved. Although 5d fields come in com- 
plete SU(5) multiplets, the orbifold projections P and P' imply that the zero modes 
under this action do not fill out full SU(5) multiplets. From a 10 we just get U and E 
zero modes, while from a 5 we get D. The remaining components of a full SM generation 
are realized at the zero mode level by taking another copy of 10 and 5 in the bulk and 
flipping the overall sign of the action of P' on these multiplets. As a result, we have 
zero modes which fill out the full matter content of T or F at the zero mode level P,^]. 
However, since the components of T or F arise from different 5d SU(5) parent multiplets, 
the interaction of X and Y gauge bosons with these 'split' multiplets do not convert SM 
quarks into SM leptons. Rather, they convert SM fermions into superheavy (m ~ 1/i?) 
exotic partner states. Only with multiple X, Y exchange can a baryon-number- violating 
operator involving only SM fields be produced. Such operators are at minimum of dimen- 
sion 8 and are irrelevant for proton decay experiments in the foreseeable future]^ Thus, 
we now focus on proton decay due to interactions of X, Y gauge bosons with SU(5)-brane 
matter. 



^If some or all of the SM fermions come from split multiplets localized in the bulk, X, F-boson- 
mediated proton decay via dimension-6 operators will occur if couplings analogous to Eqs. (^ and (||) 
mixing different bulk SU(5) multiplets are present. One can think of such terms as non-diagonal brane- 
localized kinetic terms (in a basis chosen to make the bulk kinetic term diagonal). We do not consider 
this case in detail here. 
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The only AB = 1 operators induced by single X oi Y gauge boson exchange are 

Otf = — ^Y^aihj duiU^j LiiQij (12) 

Ott = —jj^Y.^i^j'^i'^jQLiQLj- (13) 

c ij 

Here and below our Weyl spinor notation is such that ^/'Dirac = (^l, The operators 

Otf and Ctt arise from gauge exchange between T and F and between two T multiplets 
respectively. The generation indices i,j label the matter fields in the 'locality basis,' |1^ 
(in which the fields are defined by their location) so that = 1 or Oj = depending on 
whether Fi is on the SU(5) brane or not. Similarly 6j = 1, depending on whether Tj is 
localized on the SU(5) brane, or not. 

The fields in the locality basis (which is also the gauge basis) get mass via Yukawa 
interactions characterized by the matrices A„, A^, and Ag. Low-energy experiments de- 
termine only their diagonal form, A^"^^, A)^'^^, A^"^^, and the CKM matrix, Vckm^ 

Af s = LiXuRu, Af ^ = L%R,, Af ^ = L%R,, Vckm = (14) 

The diagonalizing unitary rotations of type L and R act in generation space and connect 
the gauge and mass eigenstate bases. Denoting the latter by primes, the above AB = 1 
operators take the form 

2 2 

<^TF = ^,Y.<^^h {(d^R%{^Ri)A^'LULM, - (dZRlU^Ri),{L,eML^u', 

(15) 

(note that, for our purposes, neutrinos can be treated as massless) and 

2 2 

<^tt = ^Y: hh,{^^Ri),{T^RlUKu',),{Lj,^ . (16) 

Here, in contrast to the minimal version of 4d GUTs, the coefficients of the operators 
involve combinations of elements of L and R matrices beyond those determined by Vckm- 

We start our discussion with the case considered originally in 0: all generations 
located at the SU(5) brane. We use the approximation L ~ i? ^ 1 and focus on the 
mode p vr'^e^. The relevant operator reads 



with only first generation fields. This gives rise to the decay rate (see, e.g., [pO|] ) 

Taking the hadronic parameter a = 0.015 GeV^ |2T], the pion decay constant /^r = 0.13 
GeV, the chiral perturbation theory parameters D = 0.80 and F = 0.47, the unified gauge 
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coupling gl/{An) = 1/25, and the renormalization coefficient Aji = 2.5, the resulting life 
time is ^ 

l/T(p ^ vrOe+) = 1.4 x lO^Vars x (j^i^) • (19) 



Using the Super-Kamiokande limit of 5.3 x 10^^ years |2^,|2^, this leads to the bound 
Mc > 0.8 X 10^*^ GeV § (see also Q). Given the usual unifications scale Mgut ^ 2 x 10^'' 
GeV, almost no validity range for the higher-dimensional field theory is left. 

Let us now turn to a better motivated case where only the third generation is located 
on the SU(5) brane. This scenario predicts the large mass of the third generation, b — t 
unification, and the absence of a similar mass relation in the first two generations. Given 
the observed structure of Vckm, we assume that all matrices L and R are near the 
unit matrix. Thus, proton decay will be suppressed by the small off-diagonal elements 
of these matrices. To minimize the number of such elements, is advantageous to have an 
anti-neutrino and a strange quark in the final state. The relevant operator reads 

= ^ 4^ Kid'L {R%3{RiML,hi . (20) 

The resulting decay rate, calculated using the analysis of and including the phase 
space suppression due to the mass of the K~^, is 

2\2 2 , ^9/2/i\2/ 2 



31 



2 



(21) 

with mA = 1.1 GeV and vhk = 0.5 GeV. The parametric supression by small mixing 
angles agrees with the estimate of [|lll]. If, guided by Vckm, one estimates the 2-3 and 



1-3 mixing elements by 0.05 and 0.01 respectively, one finds a life time 

l/r(p ^ K+V,) ~ 6.6 X lO^Vears x [-^o^^) ■ (22) 

Thus, to observe an effect, one requires either experimental progress beyond the currently 
considered multi-megaton detectors (see e.g. [^]), or a very low compactification scale, 
or off-diagonal elements of the unmeasured mixing matrices larger than the naive CKM- 
based estimate. 

Finally, we consider what is arguably the best-motivated scenario, namely, a model 
with only T3, F3 and F2 on the SU(5) brane. Following the discussion of |T^, one finds 



that this model provides an explanation of all fermion mass hierarchies except for the 
lightness of the first-generation up-quark and electron and the neutrinos in terms of the 
single parameter (M/M^)^/^. The dominant decay modes are p — > K'^fx^ and p — K^v^. 
For the striking K^fi~^ mode the relevant operator is 

O = (^I)l3(^«)31 . (23) 

Similarly we find that the K^v^ mode is also supressed by two 1-3 mixings.Q This gives 

^In Rcf. the modes p K^u^ and p /i+tt" were stated to be dominant with amplitudes 
supressed by one 2-3 and two 1-3 mixings and by one 1-2 and two 1-3 mixings respectively. 
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rise to the rate 



(24) 

The hfe time is 

l/T{p ^ ^ 2.1 X lO^Vears x [j^^]' ■ (25) 

Thus, in a phenomenologically well-motivated realization of this scenario, with M ~ 
10^^ GeV and [MjM^^I'^ ~ 30, the next generation of proton decay experiments could 
'discover' orbifold GUTs via the striking signal of a K^fi~^ final state. (Note that in the 
case where T2 but not F2 is on the SU(5) brane proton decay could also be discovered 
by the striking p —>■ K^^x'^ mode. In this case the mixing supression is only by two 1-2 
mixings [0.) 

We now turn to the non-minimal case where proton decay arises from the higher- 
derivative interactions of the X, Y bosons with matter on the SM brane. This is of interest 
especially in the minimal GUT models where all matter in localized on the SM brane 0. 
The proton decay rate that results from the operator Eq.(0) is 

l/r(p ^ = 3.5 X lO'W X ^ {^^)' {^JL_f . (26) 

Assuming natural 0{1) values for the unknown coefficients b and c, this leads to a 
potentially observable rate. This same mechanism can also apply to split multiplets, but 
in this case there is a further suppression of M^/M in the rate for each bulk field. 



4 Conclusions 

Proton decay is one of the most characteristic and striking signatures of conventional 4d 
GUTs. In orbifold GUTs, although dimension 5 operators are naturally absent, proton 
decay induced by the exchange of X, Y boson states is enhanced. Moreover, the mini- 
mal gauge interactions of X, Y gauge bosons only induce /S.B = AL = 1 operators for 
those matter multiplets localized to the SU(5)-invariant brane. For example, in models 
with the most realistic fiavour structure, the 1st generation is not localized to the SU(5) 
brane, and proton decay arises from the mixing of brane and bulk matter to form mass 
eigenstates. As a result, minimal X, Y gauge interactions favour unusual final states, 
e.g. p ii'^/i"'", at a potentially observable rate. Thus the rates and fiavour structure of 
proton decay processes provide important probes of orbifold GUTs. In addition to min- 
imal X, Y couplings, there can exist higher-derivative brane-localized couplings of X, Y 
gauge bosons directly to 1st generation states. Although formally suppressed by Mc/M, 
these new operators do not rely on mixing and can be the numerically dominant source 
of proton decay, with rates r(p — > vr^e"^) observable in current and future experiments. 



8 



References 



[1] H. Georgi and S. Glashow, Phys. Rev. Lett. 32 (1974) 438; 
J. Pati and A. Salam, Phys. Rev. D8 (1973) 1240. 

[2] H. Murayama and A. Pierce, Phys. Rev. D 65, 055009 (2002) 
||arXiv:hep-ph/0l08l04|| ; 

G. Aharelh, F. Ferugho and I. Masina, JHEP 0011, 040 (2000) 



arXiv:hep-ph/0007254 



R. Dermisek, A. Mali and S. Raby, Phys. Rev. D 63 (2001) 035001 



arXiv:hep-ph/0007213 



[3] 
[4] 
[5] 
[6] 

[7 



[9 
[10 

[11 
[12 

[13 

[14 
[15 



G. Bhattacharyya and K. Sridhar, |arXiv:hep-ph / 1 1 1 345 . 
Y. Kawamura, Prog. Theor. Phys. 105 (2001) 999 ||arXiv:hep-ph/0012125|| . 
G. Altarelh and F. Ferugho, Phys. Lett. B 511, 257 (2001) ||arXiv:hep-ph/0l0230l 



L. HaU and Y. Nomura, Phys. Rev. D 64 (2001) 055003 ||arXiv:hep-ph/0103125|| . 

A. Hebecker and J. March-RusseU, Nucl. Phys. B 613 (2001) 3 
||arXiv:hep-ph/010616q . 

A. Hebecker and J. March-Russell, Nucl. Phys. B 625 (2002) 128 
||arXiv:hep-ph/0107039| . 

T. Asaka, W. Buchmiiller and L. Covi, Phys. Lett. B 523, 199 (2001) 
|arXiv:hep-ph/0108021| . 



L. J. HaU, Y. Nomura, T. Okui and D. R. Smith, Phys. Rev. D 65, 035008 (2002) 
||arXiv:hep-ph/0108071|| . 

R. Contino, L. Pilo, R. Rattazzi and E. Trincherini, Nucl. Phys. B 622, 227 (2002) 



[|arXiv:hep-ph/0108102 

Y. Nomura, Phys. Rev. D 65, 085036 (2002) |^rXiv:hep-ph/010817q 

E. A. Mirabelh and M. E. Peskin, Phys. Rev. D 58 (1998) 065002 
||arXiv:hep-th/97l22l4|| . 

N. Marcus, A. Sagnotti and W. Siegel, Nucl. P hys. B 224 (1983) 159 ; 
N. Arkani-Hamed, T. Gregoire and J. Wacker, |arXiv:hep-th/0101233| ; 
D. Marti and A. Pomarol, Phys. Rev. D 64 (2001) 105025 ||arXiv:hep-th/0106256 |. 

A. Hebecker, |arXiv: hep-ph /Oil 2230| , to appear, Nucl. Phys. B. 

J. Wess and J. Bagger, Supersymmetry and Supergravity, Princeton Univ. Press, 
1983 



[16] J. Erdmenger, Z. Guralnik and I. Kirsch, |arXiv: hep-th / 0203020| . 



9 



[17] T. j. Li and W. Liao, |arXiv:hep-ph/020209q ; 

L. J. Hall and Y. Nomura, t^rXiv:hep-ph/02021071 ; 



S. Dimopoulos, D. E. Kaplan and N. Weiner, |arXiv:hep-ph/ 020213(: . 



L. Hall, J. March-Russell, T. Okui and D. R. Smith, |arXiv:hep-ph/0108161 



C. Csaki, G. D. Kribs and J. Terning, Phys. Rev. D 65, 015004 (2002) 
||arXiv:hep-ph/0107266|| . 

J. Hisano, H. Murayama and T. Yanagida, Nucl. Phys. B 402, 46 (1993) 

||arXiv:hep-ph/ 920727911 ; 

J. Hisano, |arXiv:hep-ph/0004266| . 

S. Aoki et al. [JLQCD CoUab.], Phys. Rev. D 62 (2000) 014506 
||arXiv:hep-lat/99 1102611 . 

M. Shiozawa et al. [Super-Kamiokande CoUab.], Phys. Rev. Lett. 81 (1998) 3319 
||arXiv:hep-ex/9806014|| ; 

Y. Hayato et al. [Super-Kamiokande Collab.], Phys. Rev. Lett. 83 (1999) 1529 
|arXiv:hep-ex/ 99040201 . 



Y. Suzuki et al. [TITAND Working Group Collab.], |arXiv:hep-ex/0110005| . 

R. Dermisek and A. Mali, Phys. Rev. D 65 (2002) 055002 ||arXiv:hep-ph/0108139 



10 



